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INFORMATION-Book Shelf 


D-N-LH2 


The list of new DERIVE related books is very extended this time. I want to thank many 
of you for your allusions to new publications. I specially would recommend two Ger- 
man books [7] and [8] to the non german speaking friends. I think that these two books 
are so rich in mathematical and "Derivish" contents, that you need not know the con- 
necting words in most cases. We all can be proud that all the authors - with only one 
exception - are members of the DUG. 


DERIVE-BOOK-SHELF 


[1] Linear Algebra Experiments Using the Derive Program, M. Salter & L. Gilligan, 

GILMAR PubL, P.O.Box 3676, Cincinnati, Ohio 

[2] Computer Algebra Systems in the Classroom, J.Monaghan & T.Etchells, 

Centre for Studies in Science and Mathematics Education, The University of Leeds, Leeds LS2 9 JT 

[3] Discovering Calculus with DERIVE, J. Johnson & B.Evans, John Wiley & Son, Inc. 1992 

[4] Learning Mathematics through DERIVE, l.S.Berry, E.Graham & A.l.P. Watkins Ellis Horwood. 1993 

[5] Algebraprogramme und Tabellenkalkulation im Unterricht, W.Cyrmon, Holder-Pichler- Tempsky , Wien. 1993 

[6] Mathematik mit DERIVE, W.Koepf, A.Ben-Israel u. B. Gilbert, Vieweg. 1993 (Diskette dazu erhaltlich) 

[7] DERIVE -ein mathematischer Assistent, Landesinstitut fur Schule und Weiterbildung (Kayer, Neveling) Best.Nr .1759 

beim Soester Verlagskontor, Gabelsbergerstr. 1, D-59069 Hamm (inkl. Diskette) 

[8] Mathematik mit DERIVE, Arbeitsblatter zur experimentellen Mathematik, Mauve und Moos Diimmlerbuch 4588, 

Diimmler. 1993 (vor kurzem erst erschienen. Diskette erhaltlich) 

[9] Mathematikstunden in der ll.Jahrgangsstufe, schriftl Hausarbeit von St.Ref Yasser Kanaan, 

Hans-Leinsberger-Gymnasium Landshut 

[10] Integration mit Derive, W .Koepf und A.Ben-Israel, 

Aufsatz in Didaktik der Mathematik 1, 1993. Bayerischer Schulbuchverlag 

[11] Zur Berechnung der trigonometrischen Funktionen, W.Koepf, Preprint A/1 6-93, 

FU Berlin - Fachbereich Mathematik 

[12] Ein elementarer Zugang zu Potenzreihen, W .Koepf, Preprint A/22-93, FU Berlin -Fachbereich Mathematik 

[13] Usando Derive - Fundamentos da Matematica, A. J.P. Watkins, Editora Grypho, 1993. Rio de Janeiro 


DERIVE USER GROUP UK 
Meeting 1993 


On 20th September this meeting took place in Birmingham. 38 participants followed the presentations: 
Prof. J. Berry: Derive International Journal & Conference, Dr. D. Salinger: Integrating like a 
computer algebra package. Dr .B. Kutzler: Derive in Europe, Short presentations from meeting 
attendees, K. Williamson: Finite fields alive with version 2.5, Dr. D. Stoutemyer: The Future of 
Derive and others. 

/ have to thank all members for coming. My very special thanks to Terence Etchells for organizing 
the meeting. He had a lot of work to do and to overcome a big portion of stress. At last thanks to the 
staff of the TMT Conference in Birmingham for the wonderful cooperation and the nice days in Bir- 
mingham. Thanks Pam, Mrs Vicepresident of the DUG. 
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Liebes DUG-Mitglied, 

mit groBem Vergniigen lege ieh Ihnen nun die bereits 
12. Ausgabe Ihres DNL vor. Ja, ieh sage ganz bewusst 
"Ihres" DNL, denn diese Zeitsehrift kann nur dureh Ihre 
Mithilfe, Ihr Engagement und Ihr Interesse bestehen. 
Dank Ihrer Zusehriften gelingt es immer wieder ein, 
wie ieh glaube, interessantes Magazin zusammenzu- 
stellen. 

Bitte bleiben Sie weiterhin so aktiv. In dieser Ausgabe 
konnte zum ersten Mai ein Beitrag aus Asien aufge- 
nommen werden. Afrika und Australien fehlen noeh 
auf der Autorenliste. Ieh moehte vor allem die zahlrei- 
ehen Mitglieder aus Australien ermutigen, mir zu 
sehreiben. Als Anreiz setze ieh fiir den ersten australi- 
sehen Beitrag eine Jahresmitgliedsehaft bei der DUG 
als Preis. Das gilt natiirlieh aueh fiir Afrika, wenn die 
DUG dort vertreten sein sollte. 

Ieh danke auf diesem Weg den vielen Freunden, die 
uns - meiner Frau und mir - immer wieder sehreiben. 
Ieh kann unmoglieh alle Briefe gleieh beantworten, 
aber Ihre Aufmunterungen und Anregungen maehen 
mir die vielen Stunden hinter dem Computer noeh 
erfreulieher. 

Meine Frau war der Meinung, in eine Dezemberausga- 
be gehort nun einmal ein Christbaum, so erfand ieh 
eine Christbaumfunktion (siehe naehste Seite - das 
XMAS.MTH-file ist auf der Diskette). Bei derartigen 
Problemstellungen wundert sie sieh dann, dass ieh viel 
Zeit hinter dem Computer verbringe. Ieh habe sehon 
Angst, dass ieh zu Ostem ..., na, ein Osterei ginge ja 
noeh, aber ein Osterhase??? 

Sie linden als Beilage eine Reehnung fiir Ihren Mit- 
gliedsbeitrag 1994. Ieh hoffe natiirlieh, dass Sie Ihre 
Mitgliedsehaft emeuern. Sehieken Sie bitte wie bisher 
einen Euroseheek oder geben Sie uns mit der Allonge 
Ihre Kreditkartennummer bekannt. Fiir die Mitglieder 
aus Osterreieh liegt ein Zahlsehein bei. 

Im naehsten DNL werde ieh iiber die DERIVE- Konfe- 
renz in Krems beriehten. Bei dieser Gelegenheit moeh- 
te ieh Sie auf die naehstjahrige Konferenz in Plymouth 
hinweisen. Eine Information liegt diesem DNL bei. 

Herr Dr. Wiesenbauer von der TU Wien hat verspro- 
ehen, eine Kolumne fiir Zahlentheorie und Algebra mit 
DERIVE einzuriehten. Herzliehen Dank dafiir. Wir 
sind gespannt auf seine Ideen. 

Noeh ein Gedanke meiner Frau: ieh werde im naehsten 
DNL die Gelegenheit wahmehmen, mieh personlieh 
naher vorzustellen. Ieh glaube, das gehort aueh dazu. 

So bleibt mir nur noeh, Ihnen ein sehones Weihnaehts- 
fest und ein gutes, gesundes und vor allem friedliehes 
Jahr 1994 zu wiinsehen. 


Dear DUG-Member, 

It is a great pleasure for me to present the 12th issue of 
your DNL. Yes, I insist on saying "your" DNL, beeause 
this paper ean only exist with your help, your enthusi- 
asm and your interest. With your help it is possible to 
make up an interesting journal. Please remain in your 
aetivity in the future! In this issue you will find the first 
time a eontribution from Asia. I miss Afriea and Austra- 
lia on the author’s list. I espeeially want to eneourage 
our many members from Australia to write. As an in- 
eitement I grant a one year membership for the first 
Australian and Afriean eontribution. 

By the way I have to thank the many of you, who wrote 
to me. Its not possible for me to answer all the letters 
immediately, but your friendly words and your ideas 
embellish my hours behind the PC. 

It was my wife’s idea that in a Deeembers issue a 
Christmas tree should not be missing. So I ereated a 
Xmas-tree-funetion (see next page -XMAS-MTH is on 
the diskette). She sets the problems and then she is 
wondering at the time I spend at the eomputer. I fear 
Easter time, an Easter-egg is would not be too diffieult, 
but an Easter-Bunny. ...? 

In this issue you will find an invoiee for the dues for 
1994. I hope that you will renew your membership. 
Please send a eheque or use the attaehed form to allow 
us to eharge your erediteard (VISA or Master/Euroeard 
only). 

I will give a report on the 2nd DERIVE Conferenee in 
Krems In the next DNL. At this oeeasion I want to re- 
mind you on the DERIVE Conferenee in Plymouth in 
July. You ean find a respeetive information in this DNL. 

Dr. Wiesenbauer from the Teehnieal University of 
Vienna has promised to install a speeial eolumn in the 
DNL: DERIVE - Algebra and Number Theory. Many 
thanks and we are looking forward to his future eontri- 
butions. 

Another one of my wife's ideas (she is full of them!): in 
the next DNL I will introduee myself 

At the end of this letter I wish you a merry Christmas 
and a happy, healthy and peaeeful year 1994. I hope to 
meet you all again next year . 

With my best regards as ever 



Mit den besten GriiBen Ihr 
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The DERIVE NEWSLETTER is the Bulletin of 
the DERIVE User Group. It is published at 
least three times a year with a contents of 30 
pages minimum. The goals of the DNL are to 
enable the exchange of experiences made with 
DERIVE as well as to create a group of experts 
to discuss the possibilities of new methodical 
and didactic manners in teaching Mathematics. 

Editor: Mag. Josef Bohm 
A-3042 Wtirmla 
DXust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 

Preview: (Contributions for the next issues): 

Fluid flow in DERIVE, Reuther a.o., BRA 

Applications in Electrical Engineering, Scheuermann, GER 

Minimization of a "Flat function", Lopes, FOR 

Stability od systems of ODEs, Kozubik, SLO 

Los desplazamientos en las funciones elementales, Rarnos, ESP 

Computer Algebra and Excel, Sjdstrand, SWE 

Algebraic Operations on Polynomials in DERIVE, Roanes, ESP 

Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, A 
DERIVE in Austrian Schools, some examples, Lechner & Eisler, A 
Tilgung fremderregter Schwingungen, Klingen, GER, and others 


Contributions: 

Please send all contributions to the below ad- 
dress. Non-English speakers are encouraged to 
write their contributions in English to under- 
line the internationality of the DNL. It must be 
said, though, that non-English articles will be 
warmly welcomed nonetheless. Your contribu- 
tions will be edited but not assessed. By sub- 
mitting articles the author gives his consent for 
reprinting it in the DNL. The more contribu- 
tions you will send to the Editor, the more 
lively and richer in contents the Derive-News- 
Letter will be. 




Original DOS version from 1993 versus DERIVE 6.10 version from 2006 


Impressum : 

Medieninhaber : DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrif t 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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P .Gondrie. Tilbur2. Netherlands 

... Further more a question of DERIVE: Is it possible to declare a variable of type integer. For 
example: cos (k n) becomes then (-1)'^. Now it cannot simplify this answer. 

#1 : k :€ Integer 


In 1993 it was not possible, but now you can: 


k 

#2: COS(k*n) = (-1) 


Dr.F.Schumm. Stuttgart. German- 

Beiliegend schicke ich eine kleine ausgearbeitete Erwiderung zu "Mr. Robins 2nd problem" 
(DNL#n ): 


Da Herm Robins Problem allgemeinerer Natur ist und bei DERIVE immer wieder auftaucht, gehe ich 
ausfiihrlicher darauf ein. (Because this problem is a general one I will treat it more precisely). 

1. Die "ungeschickte" Definition einer Ableitung 


Es sei: f(x) Fiir die Ableitung bilden wir DIF(f (x) ,x), also 
Geben wir nun DF ( u ) ein und wenden Simplify an, so erhalten wir wie erwartet: 3u^. 


DF(x) := — f(x) 
dx 


Das richtige Ergebnis ist jedoch triigerisch! Versuchen wir es entsprechend mit DF (2 ) , so erhalten 

wir den vielleicht nicht erwarteten Ausdruck — f (2) . 

d2 


If we define DF(x) := DIF(f(x),x), then we obtain DF(u) = 3u^ for f(x):= x^ but we find also - probably not 
d 

expected — f(2) . 
d2 


2. Begriindung fiir das "falsche" Ergebnis 

Schreiben wir DF (2 ) so meinen wir, da6 erst nach x differenziert werden und dann im Resultat x 
durch 2 ersetzt werden soil. DERIVE hingegen substituiert erst und versucht dann zu differenzieren! 
Soweit ich beurteilen kann, werden grundsatzlich bei Funktionen F(x, y, z):= bei Ersetzen der Argu- 
mente x, y, z durch andere Ausdriicke z.B. F(2, b - d) tiberall auf der rechten Seite der Funktionsde- 
finitionen sofort diese Ersetzungen durchgeftihrt, gleichgtiltig wo. 

If we write DF(2) we intend to substitute 2 forx after forming the derivative, but DERIVE first substi- 
tutes and then does the derivation wrt the "new" variable. 


3. Losung des Problems 

Wir mtissen die Reihenfolge "erst Differenzieren, dann Substituieren" ktinstlich erzwingen. Das Er- 
zwingen von Reihenfolgen lasst sich in DERIVE (fast?) immer mit dem Grenzwert DIM und einer 
Hilfsvariablen (hier x_) durchfiihren. 

Wir definieren daher die Ableitung DF(x) besser so: wir differenzieren erst f(x_) nach x_ und ersetzen 
im Resultat x_ durch x mit Hilfe des Grenztibergangs DIM x_ ^ x. 

d 

#2: DFCx) := 1-iii f(x_) 

x_-^x d x_ 

#3: [DFOO, DFC2)] = Eii , I2] 


Jetzt ergibt auch DF(2) nach Simplify auch richtig 
den Wert 12. 
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Wenn man sich die UTILITY-Files von DERIVE ansieht, so findet man an vielen Stellen 
LIM-Symbole, wo man sie zunachst nicht erwartet. Meist handelt es sich dann um das gleiche Prob- 
lem, cine bestimmte Rechenreihenfolge einzuhalten, bzw. salopp formuliert: "Erst allgemein berech- 
nen, dann substituieren." 

("You should first calculate generally and then substitute!" So we force DERIVE to do so, using UM.) 
More examples from Mr Robin 's mail - now treated correctly: 


fl: 


[/ X'DFCx) dx, J C^f + ii)'DFCx) dx] = 


4 4 

i-x i-x i 

, + II'X 

4 4 


# 4 : 


[/ X'DFCx + ii) dx, DFCx + y)] = 


3 

Cx + u) ■ Ci-x - u) 


4 


2 

3-Cx + y) 


Dr. Hubertus Wohl. Freibure. Germany 
Solar physics with DERIVE on the HP95XL 

As an astrophysicist with the actual working field of experimental solar research I am using 
since more than two years the Hewlett-Packard Palmtop PC (HP95IX) including -since sev- 
eral months - the DERIVE 2.5 program. Here are some of my experiences: 

For a mobile PC with a mass of about 300 grams, which is always with you on the road, in 
the train, and on the plane, the main usages of the HP95LX are of course the time planning 
and telephone infos. In addition many statistical reductions can be performed using the 
LOTUS 1-2-3 (vers. 2. 2) package and the included HP calculator of the standard HP95LX. 
Since the (upgraded) version has 1 MB of RAM and I use a 2 MB SRAM card in the A-drive 
there is enough storage for the DERIVE program, its additional subroutines and all user 
group programs. 

My applications of DERIVE are mainly curve fitting with determinations of the errors of the fit 
parameters, the plotting of data and curves, and applications of Fourier transforms to time 
series to determine periods. Most of these usages are during travelling, scientific confer- 
ences and observing periods at a solar observatory on the Canary Islands. In general no 
printouts are needed. Since I use at my office MAPLE V on a UNIX based workstation a part 
of the using of DERIVE is also testing of results from MAPLE applications. Of course this is 
limited to not very time consuming operations, because the processor speed of the HP95IX is 
rather slow. 

A special feature of my HP95LX is the fact that it operates with solar power: it uses two 
NiCd rechargeable batteries lasting for several hours of computing time. The NiCd cells are 
recharged with solar energy powered equipment. 


Josef Lechner, Viehdorf, Austria 

Beim Befehl ITERATES sollte der Iterationsschritt als 5. Parameter dem Benutzer zuganglich sein, 
damit wtirde der ITERATES -Befehl viel an Flexibilitat gewinnen. 

Josef would like to have an additional 5th parameter in the ITERA TES-command for the iteration step. 
ITERATES would gain much more flexibility. 

Bin Beispiel; Das Abarbeiten einer "Befehlsliste" mit verschiedenen Funktionen ware mdglich. 

As a example working through a "command list" containing various functions would become possible. 
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bliste : = [FI (x) , F2 (x) , Fn(x)] 

PROGRAM (bliste) :=ITERATES (ELEMENT (bliste, k) , x, bO , DIMENSION (bliste) , k) 

mit dem Ergebnis (with the result): 

hisher (until now): [xO, F(xO), F(F(xO)), F (F (F (xO ))),... ] 

unddarm (and then): [xO, Fl(xO), F2(Fl(xO)), F3 (F2 (FI (xO ))),... ] 

D-N-L (from 2006): ITERATES didn't change since 1993, but now we can write a little pro- 
gram which performs applies to Josef Lechner's suggestions: 


nestedCfuncs , xO, x, f) := 
Prog 
Loop 

If Din (f lines) = 0 exit 
#1: f_ := FIRST Cf lines) 

xO := LinCf_, X, xO) 
fiines := REST (f lines) 
xO 


#2: nested 


nested 


2 

X 1 


SINCx), 1 + X, 


SINCx), 1 + X, 


JT 

2 } 


, t 


SINCt) + 2-SINCt) + 1 


The argument funcs is the list of functions which should be nested, followed by the initial 
value xO for variable x (by default). Compare the two next function calls: 

3 

#4: r 2 x-yl y -COSCy) 

nested(LcOSCy) 1 x-y , e J , ff) = e 

2 

#5: r 2 x-yl x 

nested(LcOSCy) 1 x-y,e y) = e 


Klaus Herdt, Osnabrtick, Germany 

Wie bereits in der Beantwortung Ihres Fragebogens fehlen m.E. manchmal Hinweise mehr techni- 
scher Art. Beispielsweise konnte ich mir vorstellen, dass es eventuell auch anderen Anwendern nicht 
immer gelingt, die auf der Hardware ihres Rechners prinzipiell verfiigbaren (und unter anderen Pro- 
grammen auch durchaus einsetzbaren) Grafikmoglichkeiten auch unter DERIVE zu nutzen. Ich habe 
beispielsweise auf meinem PC cine Grafikkarte EIZO VA4I (I MB) im Einsatz, die theoretisch cine 
Auflosung von 1024x768 bei 256 Farben ermoglicht. Die vom Hersteller angegebenen Modi 37h, 
bzw. 38h (16, bzw. 256 Farben bei einer Auflosung von 1024x768) werden von DERIVE zwar syn- 
taktisch akzeptiert, nicht jedoch unterstiitzt. Neben DERIVE benutze ich bisweilen auch MAPLE V, 
rev.2 und dieses System bedient sich bei Betrieb unter DOS ganz analog zu DERIVE einer Initialisie- 
rungsdatei (MAPLE.INI). Hierin lassen sich die Darstellungsarten 800x600x16 bis 1024x768x256 
ansprechen durch die (genormten?) Modusangaben 258 bis 261 (die dann wohl intern in die herstel- 
lerspezifischen Kennzahlen umgesetzt werden). Naturlich liegt es nahe, auch bei DERIVE einen Ver- 
such mit diesen Modi zu unternehmen. Line Anfrage beim Hersteller erbrachte lediglich die wenig 
befriedigende Antwort, dass dann wohl die Hardware nicht korrekt arbeite, was durch das MAPLE- 
Beispiel ja wohl widerlegt sein diirfte. 

Vielleicht hat der eine oder andere Anwender ahnliche Erfahrungen gemacht und das Problem losen 
konnen?! 


p 


6 
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D-N-L; I hope that there are any members who are willing to offer their experiences with "trouble 
shooting" graphic adapters or with trouble shooting DERIVE in connection with graphic cards. I have 
similar problems not with DERIVE but with GRASP and a Cirrus Logic Card. 

Wolfgang Propper, Nlirnberg, Germany 

Auf beiliegender Diskette befindet sich ein Icon, das ich mir letzte Woche gestrickt babe. Ich stelle es 
geme der Allgemeinheit zur Verfiigung, d.h., ich verzichte auf alle, mir aus der Erstellung dieses gros- 
sen Werks zustehenden Urheberrechte. Damit mtissen DERIVE-Anwender, die von WINDOWS aus 
starten, nicht mehr auf ein eigenes Icon verzichten und konnen DERIVE wie andere schone Pro- 
gramme (MathCad, Mathematica} durch Doppelklick starten. 

Die Installierung ist denkbar einfach: DERIVEXM.ICO wird ins DERIVE- bzw. ins WINDOWS-Ver- 
zeichnis kopiert. Im ProgMan wird im Mentipunkt Datei-Eigenschaften . . . der Pfad fur DERIVE und 
gegebenenfalls das Arbeitsverzeichnis angelegt und dann mitAnderes Symbol. ... \md Durchsuchen 
dieses Icon mit DERIVE verbunden. Mitgliedern der DUG konnen Sie diese Datei auf der 93er 
DNL-Datei mitschicken. 

Noch viel schoner ware es allerdings, wenn moglichst bald eine WINDOWS- und/oder eine OS/2- 
Version von DERIVE erschiene. Am Icon muss es nicht mehr scheitern. 

D-N-L: Many thanks Mr Propper for this nice Christmas gift for the User Group. With pleasure I 
copy your creation DERIVEXM.ICO onto the disk of the year 1993. But that our users who don't 
work with the XM-version have not to miss their icon I tried - stimulated by your idea - to "paint" 
DERI. ICO and DER2.ICO. Installation is very simple. Copy the icons into the DERIVE- or WIN- 
DOWS -directory and connect DERIVE with one of the icons using the Program Manager. If anyone 
of you have created other DERIVE-Icons, then please let me know. Let's have a DERIVE-Icon-Exhi- 
bition. Merry Christmas, Wolfgang, and many Christmas presents for you. 



Hubertus WohL Freiburg, Germany 

Zu den Grafiken Erkki Ahonens sollte deutlich darauf hingewiesen werden, dass die Resultate im We- 
sentlichen auf Grund eines Programmfehlers entstehen und bei vielen DERIVE-Nutzem (so auch bei 
meinem Rechner mit Version 2.5) gar nicht verifiziert werden konnen. 

(Mr Wohl wants to point to the fact that the results of Erkki Ahonen are caused by a bug in 
the program and that many users who are working with later versions cannot verify them, 
ed.) 

Gunter Scheu, PfinztaL Germany 

Hier sind einige Befehle zum Halbierungsverfahren, wie wir es im Unterricht durchgeftihrt haben. Die 
Anregung erhielten wir im DNL#1 1 (Dyer). 

(Some functions to the Bisection Method as we used them in school. We were inspired by Mr 
Dyer's submission in DNL#1 1 .) 
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#1: Das Ha1 bier lings verfahren - Bisectioniiethod - G .Scheu ,HALB_A1 .HTH 

# 2 : 


a 

II 

b 


a + b 


a 

2 

b 

fCa) 

fCii) 

fCb) 

Till f 

Till f 

Till f 

x-ha 

x-^Ca + b)/2 

x-^b 


: HCf, a, b) := 


#4: Sclirittweise Berechnung - Step by Step Calculation 


#5: HCx - 4, 1, 4) 


# 6 : 


#7: 


# 8 : 


4 


a 

n 

L 

1 

1 

5 

A 

[ 

fCa) 

2 

fCn) 

fCb) 

-5 

9 

12 

' 2 

4 

4, 1, - 

5 " 


a 

II 

2 . 

b 

1 

7 


5 

fCa) 

4 

fCn) 

2 

fCb) 

-5 

15 


9 


16 


4 


# 10 : 


17 


#11: H 


fCa) fCn) fCb) 

15 ll 9 

16 64 4 

2 


# 12 : 


a 

1.75 

fCa) 


-■-] 
4 fi J 

II 

1.9575 

fCii) 


b 

2.125 

fCb) 


-0.9575 -0.24609575 0.515625 


[j^^HProqraiiiiierte Berechiiiinq - Calculation proqraiiiied 


HiCfi a I b, eps) := 

If LIHCfi X, + b)/2) = 0 V ABSC^ - b) < eps n 
HCf, a, b) ^ 

#14: If LinCf, X, a).LinCf, X, Ca + b)/2) < 0 

HiCf, a, Ca + b)/2, eps) 

HiCf, Ca + b)/2, b, eps) 


2 

#15: HiCx - 4, 1, 3, 0.001) 


a n b 


#16: 


12 3 

fCa) fCni) fCb) 


-3 


0 


5 


p 
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2 

#17: HiCx - 4, 1, 5, 0.001) 


a II b ' 


#18: 


12 1 
fCa) fCni) fCb) 


I -I 0 5 J 


2 

#19: HiCx - 4, 1, 4, 0.001) 


# 20 : 


1.999755859 

fCa) 


9.765625-10 


II b 

2 . 000122070 2 . 000488281 

fCii) fCb) 

-4 16 

4.8828125-10 

8191 


Hellmut Scheuermann, Hofheim/T., Germany 

We have found a bug? In version 2.52, when we wanted to show that abs(x) cannot be differentiated at 
xo = 0, what to do with the result SIGN(O)? 

#1: FCx) := |x| 

FCx + li) - FCx) 

#2: RSCx) := 

h 

#i: L"ini"it fro™ the left: 

BTLCx) := 1"ini nS(x) 

#4: IwO- 

#5: PITLCO) 

# 6 : -1 

#7: wie erwartet rnclit-ig - how expected correct! 

HTRCx) := 1"iii flSCx) 

#8: IwOe 

#9: riTRCO) 

# 10 : 1 

#11: auch r-ichtng, aber jetzt - correct, too, but now 

riTBCx) := 1niM HSCx) 

#12: IwO 

ill: HTBCO) 


IGN(O) 


The recent version of DERIVE gives an answer for sign(0): 

#16: PITBCO) = ±1 
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Message 3004 was entered on 5124193 at 8:52 AM. 

From ROGER FOLSOM to SOFT WAREHOUSE about #2974 / COMMENTS ON DERIVE 

I’ve discovered that one can get around the 512 limit on author line length by using the #number facility. That’s 
worth mentioning in the manual. 

Message 3019 was entered on 5/28/93 at 2:10 AM. (Read 10 times) 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 

BUILD. A brief paragraph listing the advantages of BUILD over Author (bypassing line edit buffer size of 512 
bytes, recording the sources of an expression, etc.) would be well worth including in the manual. 


Message 3020 was entered on 5/28/93 at 2:11 AM. (Read 10 times) 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 

POLYNOMIAL ORDERING 

You write” I think we all agree that the polynomial a + b X + c should be displayed in the above order. ”I for 
one most defmitely do not agree. In at least some areas of applied mathematics (e.g. my own field of economics) 
the usual is to display polynomials in ascending powers, e.g. as ao + ai x + a 2 x^ and so forth, simply because the 
degree is often not known initially, and adding or subtracting a term is a notational pain because it requires rela- 
beling the coefficients. Turning your polynomial into a cubic would make it 

a x^ + b x^ + c X + d or ao x^ + ai x^ + a 2 X + as; 

note in your version that coefficient a no longer is for x^, and in the last version that the coefficient subscripts 
do not match the term exponents. My apologies for belaboring the obvious, but such issues are what drives 
some of us to prefer ascending powers. [If the model will be estimated statistically, the degrees of polynomials 
used may depend on the data set and the degrees of freedom required — more data, or ’’tighter” data (around a 
regression hypersurface), lets one consider higher degree polynomials; a model may go through several itera- 
tions in which polynomial degrees are varied depending on the statistical results]. Admittedly, ascending pow- 
ers makes division of polynomials awkward (which is probably why ’’real” mathematicians such as yourself 
prefer descending powers), and in a division situation even I would order a polynomial in descending powers 
(as you would), but for many of us that case is an exception (the last time I divided any polynomials was in the 
early 1970s). 

You write further: ’’Naturally, a switch could be added that would reverse the order of all sums, but 
this radical approach would screw up the display of quadratics like the above.” For me, a quadratic or cubic 
displayed as 

a + b x^ + c x^ + d X or ao x^ + ai x^ + a 2 x + a 3 

is NOT ’’screwed up” — it is displayed normally. That simple switch is what I, and I think others, want. 

I won’t think anyone is asking DERIVE to remember whether a polynomial came from a Taylor ex- 
pansion. We just want the option to force Derive to order all polynomials in ascending powers, with the con- 
stant term listed first. If you want to get fancy, you might include an override to cause polynomials in a par- 
ticular expression (or subexpression) to be ordered opposite to however the global ’’polynomial ordering” 
switch is set. 

AUTO-SIMPLIFICATION. 

Thanks much for inviting our attention to the undocumented feature (entering a Transfer Load or Transfer 
Merge filename with a Ctrl-Enter, or Ctrl-J) that causes every line to be simplified as it is read in. I have dis- 
covered this to be helpful, but often one wants a line entered and then simplified, to keep a record of the steps 
in a calculation. Of course one can enter each command twice: once surrounded by quotes so that Derive takes 
it as a comment, and second as a real command that Derive will simplify automatically (using the Ctrl-Enter 
option). To have a simplify command that simplified the preceding expression would be somewhat more con- 
venient. Nevertheless, the Ctrl-Enter option will be useful; thanks for un-undocumenting (aka documenting) it. 
My wife says I should ask: what other undocumented goodies are buried in Derive’s depths? 
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Message 3021 was entered on 5/28/93 at 2:13 AM. 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 
EXPAND and FACTOR 

I realize that Expand can expand a sum to be not over a common denominator. However, in doing so it often 
rearranges variable ordering (of sums and products) in ways not to my liking. I realize also that this is sup- 
posed to be fixable by declaring secondary and tertiary etc. expand variables, but I have yet to figure out how 
those work. The manual (2.08) needs work here. Some way to factor an expression out of another expression 
would be very useful; e.g. if the expression a+b appears in several terms of a sum, it would be nice to be able 
to factor it out. Also, sometimes it improves readability to be able to factor an expression out to the front, or to 
the rear; to be able to choose would be a convenience. 

Finally, for me the FACTOR command does not always work: I had a variable dyda (word mode) that I wanted 
to pull out of a two-term expression (where dyda appeared in each term), but Derive insisted on leaving it 
where it was. That is, factoring dyda out of dyda^X\ + dyda^X 2 , where X\ and X 2 were two expressions, re- 
mained as dyda^Xi + dyda^X 2 , rather than as dyda^{X\ + X 2 ). The error here probably was mine, but I donT 
know what my error was. 

Message 3022 was entered on 5/28/93 at 6:35 AM. (Read 11 times) 

From JERRY GLYNN to PUBLIC about ORDER AND EXPAND 

Inspired by R Folsom’s questions: If I aske Derive to EXPAND (a+b)^ could I get a^ +2ab +b^ and if I EXPAND 
(b+a)^ get b^ + 2ba + a^? 

Message 3032 was entered on 6/2/93 at 10:27 PM. 

From SOFT WAREHOUSE to ROGER FOLSOM about #3019 / DERIVE SUGGESTIONS 

Thanks for your wide ranging thoughtful and insightful suggestions for improving DERIVE. The following is in 

response to some of them: 

SCREEN PRINTING: When printing screen images, DERIVE actually reads the screen pixels to determine 
what to send to the printer. It is impossible to read what is not actually on the screen. However, some super 
VGA graphics adapter cards have video modes that can display 132 character columns, instead of the normal 
80 columns. If then such a mode is properly supported by the BIOS and you know the number of the mode, 
select "Other" in the Adapter field of the Options Display command and enter the video mode number. Then, 
you can print screen images containing much longer expressions. 

DOMAIN DECLARATION SYNTAX: Yes, the Declare Variable Domain command will create an expression 
declaring a variable’ s domain, just like the Declare Variable Value command creates an expression declaring a 
variable’s value. 

VARIABLE ORDERING: Using an external text editor, you can make one line INI files that just contain a 
variable ordering list. Then use the Transfer Load State command to load these truncated INI files. It should 
work. 

MATCHING PARENTHESES: Being a long-time LISPer I can certainly appreciate the advantage of jumping 
to matching parentheses. Also blinking matching parentheses would be nice. 

BUILD COMMAND: In the next printing of the manual, I will discuss the benefits and drawbacks of the Build 
command. Aloha, A1 Rich, Soft Warehouse, Inc. 

Message 3033 was entered on 6/2/93 at 10:28 PM. 

From SOFT WAREHOUSE to ROGER FOLSOM about #3020 / DERIVE SUGGESTIONS 
POLYNOMIAL ORDERING: Your points are well taken. I will look into the ramifications of a flag for revers- 
ing sums. 

VARIABLE ORDERING: Let’s see if the above change is sufficient before making additional changes. 
AUTO-SIMPLIFICATION: The Transfer Demo command reads in and displays an expression, and then sim- 
plifies it. However, it pauses after each simplification. The obvious solution is to have an option to not pause. 
Aloha, A1 Rich, Soft Warehouse, Inc. 
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Message 3034 was entered on 6/2/93 at 10:29 PM. 

From SOFT WAREHOUSE to ROGER FOLSOM about #3021 / DERIVE SUGGESTIONS 

EXPAND and FACTOR: What is the specific example where dyda will not trivially factor out of 

dyda^X\ + dyda^X 2 l 

Message 3037 was entered on 615193 at 5:56 AM. 

From JERRY GL YNN to SOFf W AREHOUSE about #3032 / DERIVE SUGGESTIONS 

Hooray for A1 Rich., .just when I think IVe wasted my breath some months ago I learn that I didn’t. This will 
inspire new requests! 


Message 3045 was entered on 6/6/93 at 4:41 PM. (Read 17 times) 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 
A1 Rich: a few responses to your much appreciated responses to my suggestions: 

VARIABLE ORDERING: you note that to get a variable ordering list appropriate for a particular paper one is 
working on, one can load an *.ini file that contains only a one line ordering list. Thanks much for the sugges- 
tion. The manual might usefully mention that loading a second *.ini file does not dump all defaults previously 
loaded by an earlier *.ini file or while within Derive; only the new file’s explicit changes are loaded. 

Incidentally, my experience (Derive 2.08) is that the ordering list cannot be more than 128 characters; I bumped 
this limit the other day with a bunch of multi-character variables. Derive’s error message was something akin to 
’’can’t fmd the end of file; contact Soft- Warehouse, ” but pruning the ordering list below 128 characters solved 
the problem. A 256 or even 512 character limit would be useful, but I realize that there are severe memory 
constraints and tradeoffs (at least in ’’ordinary” rather than XM Derive). 

I will continue to observe carefully how Derive responds (in sums and products) to my ordering list, and I hope 
I can do better at predicting and forcing the patterns I want. You are right that an ascending powers switch will 
do much to make me even happier than I now am with Derive. 

Nevertheless, I tentatively still believe that an order list (maybe call it a sequence list) going from least to most 
main, together with a ’’Derive, don’t get too clever” switch, telling Derive to order all factors (and maybe equal 
power sums) in that sequence regardless of what Derive thinks is a parameter and what is a variable, would be a 
useful option. 

EXPAND AND FACTOR: I should have explained more fully. The expression dyda^Xl + dyda^^Xl initially 
involved a common denominator, and was written as 

dyda(X\ + N^X2) 

N 

I didn’t want it over a common denominator, so I expanded it around V, and got 


2 ) 


dyda * XI 
N 


+ dyda * X 2 


Then when I simplitied (2), I got (1) back again. When I factored dyda out of (2), I got (1) back again. So I 
gave up! I wanted, but was unable to get. 


3) 


dyda 



As noted earlier, also useful would be the ability to factor an expression out of another expression; perhaps 
when the factoring command asks for the variable to factor, an expression (rather than a simple variable) could 
be typed, or selected by highlighting. 

AUTO-SIMPLIFICATION. The ’’control-enter” option for loading a file, so that the file simplifies as it loads 
is useful, but it does not substitute fully for a statement within the file to simplify a particular statement. Some- 
times one does not want all statements simplified. 

Example: when I used ’’control-enter” (or Load Demo) to load a file containing the definitions y(Nt,LNt,t,a): = 
and Y : =y(Nt,LNt,t,a), Derive simplifies these to y(Nt,LNt,t,a) and y(Nt,LNt,t,a), and then doesn’t know what 
to do with later expressions such as DIF(Y,x) because it doesn’t know what Y is. I don’t see how to use Transfer 
Demo (with or without a pause, especially since the pause comes after the simplification) to deal with this 
problem. Thus a ’’simplify previous (or next) line” command still appears useful. With it. Derive would add 
some of the advantages of a ’’compiled” program, that one could ’’run” on a slow machine while one does some- 
thing else (e.g. eat dinner, or go to bed). But simplifying everything is too much. 
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Message 3056 was entered on 6/12/93 at 5:52 AM. 

From HARALD LANG to JEFF I about #3049 / TO: Y’ ALL 


In addition to Jerry Glynn ’ s suggestion, you need to tell the printer to use the IBM character set. Furthermore, 
the output gets much nicer, in my opinion at least, if you choose a line spacing of 12 lines per inch (or maybe 
even 16). I send the following sequence of control codes when I print to laser printers from Derive: 

< ESC > E < ESC > (lOU < ESC > &1 12D 

In addition, I send control sequences to set margins and so on, but since my European paper format is different 
from yours, I suppose you will find them useless. Note: < ESC > stands for the ASCII character 27. 
Regards -Harald Lang 


Ebene Algebraische und 

Transzendente Kurven (2) 

Thomas Weth, Wurzburg, Germany 

Die Strophoide 

In einem Briefwechsel zwischen F. de Verdus und Torricelli vom 15. Marz 1645 und dem 19. Mai 1645 wird 
cine "Linie" beschrieben. die" ... in Frankreich Ala oder Pteroide genannt wird ...". Die Erfindung der 
Strophoide, wie die Kurve heute allgemein genannt wird, geht somit auf einen franzosischen Mathematiker, 
wahrscheinlich Roberval zuriick. Die urspriingliche etwas komplizierte raumliche Konstruktion geschieht mit 
Hilfe von Brennpunkten von Kegelschnitten. In der mathematischen Literatur fand die Strophoide dann 1757 
im 

IV. Band der Instituti Bononiensis Comentarii in den zwei Denkschriften von Gregor Casali De conicarum 
sectionum focis (Uber die Brennpunkte von Kegelschnitten) Beachtung, geriet aber zunachst in Vergessenheit. 
Casali schlug den Namen Pteroides torricellanea vor, um die Verdienste Torricellis (einem Schuler Galileis) 
um die Kurve zu wiirdigen. Ab Mitte des 18. Jhdts beschaftigten sich u.a. de Moivre, der eine Konstruktion in 
der Ebene fand und die Quadratur der Kurve angab, Maria Gaetana Agnesi (1748), Cramer (1750), Ricati und 
Saladini (1765) mit dieser Kurve. 1819 wurden die Kurven von A. Quetelet in dQX Dissertatio de quibusdam 
locis geometricis (Abhandlung iiber einige geometrische Orter) als Fokale bezeichnet. Durch eine Preisaufga- 
be aus dem Jahre 1840, die den franzosischen Lyceen und 1861 im Examen far die Zulassung zur Ecole nor- 
male in Paris vorgelegt wurde, entstanden zahlreiche Untersuchungen der Kurve. Zum erstenmal tritt der Na- 
me Strophoide in den Nouv. Ann. V., 1846 in der Arbeit La strophoide von Montucci auf; abgeleitet ist er vom 
griechischen axpocpo^ (= Band, Seil). Bis zum Ende des 19. Jhdts trat die Kurve unter den Namen Kukumaeide 
(Bezeichnung durch Uhlhorn - 1809 - und Lehmus - 1842), Logocyklika (Booth - 1858) und harmonische 
Kurve (Rummer - 1868) auf. 

Unter den vielen interessanten Eigenschaften der Strophoide seien hier genannt: Die Strophoide ergibt sich als 
FuBpunktkurve einer Parabel, als Bild einer gleichseitigen Hyperbel bei der Inversion am Kreis oder als sog. 
"Gegenkurve" einer Geraden. AuBerdem ist sie eine anallagmatische Kurve, also eine Kurve, die unter der 
Inversion am Kreis auf sich selbst abgebildet wird. 

The Strophoid 

In a correspondence between F. de Verdus and Torricelli (1645) a "line" is described which " ... in 
France is called Ala or Pteroide ...". The strophoid was discovered out in France, probably by Rober- 
val. The first construction used the foci of conics in space. Later we can find the curve in two papers of 
Casali. He proposed the name Pteroides torricellanea. In the 18th century Moivre, M. G. Agnesi, 
Cramer, Ricati and Saladini dealt with this curve. In the 19th century especially french mathematicians 
investigated the curve. In Montucci's publication La strophoide from 1846 you can find the curve's 
name the first time. The name is derived from the greek word axpocpo^ (= band, rope). Till the end of 
the 19th century there were some other names used as Kukumaeide, Logocyklika and harmonic 
Curve. 

Among the many interesting properties of the strophoid you can find: it is pedal curve of a parabola, 
image of an equilateral hyperbola using the inversion on a circle or the so called "anticurve" of a line. 
The strophoid is an anallagmatic curve, i.e. it is mapped onto itself using the inversion on a circle 
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Die Konstruktion der Kurve 

Gegeben ist ein Kreis k(0,r), sowie zwei senk- 
rechte Durchmesser SOX and 1. Durch S 
zeichnet man einen Strahl, der 1 in A und den 
Kreis in B schneidet. Auf dem Strahl tragt man 
von S aus die Strecke AB ab und erhalt den 
Endpunkt P. Dreht man nun den Strahl um S, 
so beschreibt P eine Strophoide. 


Construction of the Curve 

Given is the circle k(0,r) and two orthogo- 
nal diameters SOX and I. Draw a ray 
through S giving the intersections A with I 
and B with k. Using the distance AB = SP 
you will obtain point P. Rotating the ray 
round S point P will describe a strophoid. 


Compare the sketch from 1993 
with the pretty GeoGebra con- 
struction from 2006. 



m GeoGebra - strophoid. ggb 


File Edit View Options Window Help 



Herleitung der Kurvengleichung (Derivation of the equation of the curve) 

Gesucht sind die Koordinaten des Strophoidenpunktes P Aus dem Strahlensatz folgt 

AO : SO = y : X, daher AO = Da nach der Konstruktion SP = AB, sind die Dreiecke ASHP 


und AACB kongruent. Damit gilt auch SH (= x) = AC = OK. 


Fur die Koordinaten von B ergibt sich somit B 


f y.y ^ 

x^r ,- — + 7 


Da B ein Kreispunkt ist, gilt: 


(x + r-r)^ + 


y -r 
— + 

V X j 


= . Fur den Strophoidenpunkt gilt dann (nach Umformung mit 


DERIVE): 


j^^(x + r) + x^(x-r) = 0 bzw. y = ±x. - — 

Vr + x 


Damit sind die Strophoiden Kurven dritter Ordnung. (The strophoids are algebraic curves of 
order 3.) 


Mit X = p cos (p und = p sin cp erfolgt mit einiger Umformung {DERIVE) die Polardarstellung der 
Strophoide und eine mogliche Parameterdarstellung (Polar form and one possible Paramter 
form): 


cos2(d , 

p = r bzw. 


x = r— 

t^+l 

y = rt— 

t^+lJ 


cos (p 


(See also the following DERIVE file.) 
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Kurzreferenz fur DERIVE User: 

Die Strophoide / The Strophoid 



Polarform: 

Parameterform: 

Algebraische Kurvengleichung: 


p = r-C0S(2-<p)/C0S((|)) 

[r.(t^2 - l)/(t^2 + 1), r.t.(t^2 - l)/(t^2 + 1)] 
y^2 (x + r) + x^2-(x - r) = 0 


The area of the loop is 2r^ — the area between the curve and the asymptote x = -r is 2r^ + 
#1: CaseRode := Sensitive 


r K 


#2 

#4 

#5 

#6 

#7: 


2 2 2 
k := C^f - f’) + y = 

1 := X = r 

ray := y = TAN(f).x 
A := [r, TAN(f).r] 

SOLUTIONS Ck A ray, [x, y]) 
0 0 


2 . r . COS Cf ) 2 . r . S IN Cf ) ■ COS Cf ) 

#«: B := [2.p.C0SCf) , 2 ■ r -SINCf) -COSCf)] 

#9: d := |A - B| 

r 


d := 


2-r.COSCf) 


# 10 : 


COSCf) 
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# 11 : 


# 12 : 


fll: 

#14: 


2 2 2 

SOLUTIONS Cray ax + y = d , [x, y]) 

2 

2.r.C0SCf) - r 2 ■ r .SIN(f) -COSCf) - r-TANCO 
2 

. r - 2- r -COSCf) r-TANCO - 2 - r -SINCO 'COSCO . 


P := L2-r-C0SCO - r, 2 - r -SINCO 'COSCf) - r-TANCO 


2-r-COSCt) - r, 2 - r -SINCt) -COSCt) - r-TANCt) 


#15: 


|^table|^[: 


2- r -COSCt) 


r, 2-r 


SINCt)-COSCt) - r-TANCt)], t, 0, n, 3] 

400 )) 



The strophoid on the Voyage 200 



It is important to work with the vector AB to 
translate the origin on the ray OB for obtain- 
ing point P. Then we can either trace point P 
or generate the locus. (I didn't label the 
points on the Tl-screen) 


#16: L2t.COSC 4>) - r, 2.r.SIN(<p) -COSClf) - r-TANCqi) 


#17: 

#18; 

#19: 

#20; 

# 21 : 


Trngononetry := Expand 

2 r-SINCip) ■ 

2.r.COSC4>) - r, 2.r.SIN(4») .C0SC4>) 

L C0SC(|>) J 

Trigononetry Auto 

r 2 r-SIN((|l) 

SOLUTIONS p. COS (4i) = 2.r-C0SCl|>) - r a p-SIN((|») = 2-r.SIN((|l) -COSClIi) 

L COSCip) 


[[ 


2-r-C0SC<p) 


— ■ "11 

cosc<p) JJ 


[p, x] 


This is the polar form; plot it changing the system of coordinates 


#22: 2-r-COSC<p) 

cosc<*>) 

Bring over one common denominator and then Trig Collect 

2 

2-r-COSCiIt) - r 

#25: 

COSC4t) 

#24: Trigonometry := Collect 

r-C0SC2-(|») 

#25: 

COSC^t) 


#26: 


2 

X = 2-r-C0SC*I») - r, y = 2 - r -SINC<^) 'C0SC<tt) 


r-SINC^f) 
COSCf) - 


pi 6 
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#27: 


#2ft: 


2 r.SINCACOSCt)) 

* = 2.r.COS(«:OSCt)) - r, y = 2 ■ r .SIN(AC0S(t)) .COSCACOS(t)) 

C0S(AC0S(t)) 


2 2 
X = r.(2-t - 1), y = r-40. - t ) 


f-4-l 


#29: SOLUTIONS Cx = r-C2't - 1) , t) 


#30: 


42.^ 


I 


r 

L r J 

L . J 


2 


2 




5 


This is the plot of the strophoid using a tool from from 
DNL#63 (2006). 




2 

X - (4 - x) 

, - 4 , 4 , - 4 , 4 

X + 4 


The area included by the bend and the area included by the 
curve and the asymptote x = -r can easily be computed by 
the CAS. 



-4 
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The red coloured area: 


The blue coloured area 


r 

f 


) 

dx 

1 r . 

1 ^ 

0 



. n J 


2 

r ■ C4 - it) 


2 




2- 


0 

y^-4 


) 

Hx 

1 ^ 

-r 




2 

r ■ Cff + 4) 


2 



As a appendix (2006) I'll show how to create the strophoid as pedal curve of a parabola. The 
pedal point is the origin of the directrix (=P). To make calculation easier we set the parabola 
into the plane that the y-axis is the directrix, i.e. the origin is the pedal point: = 4r (x- r). 

The pedal curve is the locus of all intersection points between the tangents of the parabola 
and the perpendicular lines through the origin wrt the tangents. 

This is a nice - and not too difficult - task for students. Introducing slider bars adds some 
animation. 


# 1 : 

t2-. 

fi\ 

#4: 

#5: 

fG: 


2 

y = 4 ■ r ■ Cx - r) 

[ 0 , 0 ] 

r 2 

tang := y - t) + 2-4(.r-t - r ) 

2 

^Cr-t - r ) 
x.^r.Ct - r)) 

perp :z y = - 

r 

(SOLUTIONS Ctang a perp, [x, y])) 

1 

r.(2.r - t) (t - 2.r).^Cr.(t - r)) 
t t 


Insert slider bars for r and t and see the tangents, 
perpendicular lines and their intersection points 
moving around. 

#6 is - one possible - parameter form of the 
strophoid. For plotting purposes one has to sub- 
stitute for t, because t is used as slider bar vari- 
able. Elimination of the parameter leads to the 
implicit form of the strophoid. 


#7: 


r.(2.r - u) (u - 2.r).-/Cr.Cii - r)) 

X = , y = 


#a: SOLUTIONS x = 


r r.(2.r - u) Y\ 2-r 

r = — ^ — ■ “JJ = T7- 


(II - 2 .r). 4 /Cr.Cii - r)) 2-r 

#9: SUBST y = , u, 


#10: y = - x.SIGN(r) 


#11: y = - x.SIGN(r) 


4 ^] 


2 X ■ (r - x) 
#12: y = 



Explore Chaos with 
NEWTON and RAPHSON 


Josef Bohm, Wtirmla, Austria 

I am sure that in a very near future many classical items of school mathematics will be obsolete 
if we are able and willing to use Computer Algebra tools. They will not be obsolete at all, but the way 
we teach and use them will have to change. Lets have a glance to some of teachers' most beloved 
“children”: discussion of a function - calculating the zeros, differentiating and solving equations to 
find maximum and minimum values and inflection points and using the results for sketching the 
graph. Now a simple keypress on two buttons gives the graph within less than a second and DIF in 
connection with SOLVE will do the work for most of the examples which can be found now in our 
students' textbooks. I really don t believe that we will not have to teach this any longer, but there are 
new and fascinating possibilities and we can give traditional problems some new qualities. I remem- 
ber with great pleasure Paul Drijvers' lecture held in Krems this year. He showed examples how the 
investigation of a function's properties can be changed in a thrilling way for students - and teachers 
as well. 

Another problem has a fixed place in our school mathematics. We are glad to demonstrate the 
Newton-Raphson algorithm as a first application of calculus. But now, which need is for this algo- 
rithm any longer, when pressing L in approximate mode solves the equation (DOS-DERIVE). We 
should use the occasion to focus upon the algorithm and specially on the importance of the initial 
value or initial guess of this iterative process in order to understand some “strange” results given by 
CAS-programs. 

Using DERIVE ’s graphic features it is easy to create a function for plotting the family of tan- 
gents emerging from any given initial value. So you can observe the iteration and convergence proc- 
ess. Zooming into the neighbourhood of the zero gives more insight. 

We will have a very special view on this method: 




First of all load the utility file RAPHSON.MTH, edit the function u : =f (x) and plot 
TAN_FAM(u , xO , n), which gives a family of n tangents (linear approximations) with initial value xO. 
Set Points Connected (and Small Size) and you will see the first n steps of the linear approximation 
process - with hopefully visualization of the convergence towards a zero of the function. 


#1: 

NEW_R0W(u, X, xO, n) := ITERATES(x - u/0(u, x) , x, xO, n) 


#2: 

NEW_C0L(u, X, xO, n) := NEW_R0W(u, x, xO, n) ’ 


#3: 

NEWTONCu, X, xO, n) := ITERATE(x - u/0(u, x) , x, xO, n) 


#4: 

LINE(x0, yO, dl, x) := dl*x + yO - dl*x0 


#5: 

TANGENTCu, x, xO) := LINE(x0, LIM(u, x, xO) , LIM(0(u, x) , 

X, xO), x) 

#6: 

VERTCu, xO) := [xO, 0; xO, LIM(u, x, xO)] 


#7: 

TAN_FAM(u, xO, n) := [VECTOR([VERT(u , (NEW_ROW(u, x, xO, 
X, (NEW_R0W(u, X, xO, n)) 4 ,k)], k, n)] 

n))ik), TANGENTCu, 

#8: 

NEW_ATTR(u, start, end, inc) := VECTOR(NEWTON(u , x, start 
end - start, inc) 

+ k) , k, 0, 

#9: 

NEW_ATTR_GR(u , start, end, inc, n) := VECT0R([ start + k, 
start + k, n)], k, 0, end - start, inc) 

NEWTONCu, X, 


Functions #1 - #7 are from an excellent new DERIVE publication (Book Shelf [7]). Expression #2 
returns the values of the algorithm as a column vector. 

Example 1: Take the simple cubic equation - x = 0. What can be done? It is clear that there are 
three zeros: -1,0 and +1. There is no need to apply a numeric process to find the zeros. But it might be 
interesting for students that three initial values which are very close like xq = 0.44720, xq’= 0.44725 
and xo = 0.44730 will converge towards quite different zeros. Furthermore they will recognize that 
they will need more than 10 iteration steps in case of the third initial value for noticing any conver- 
gence. 15 steps ore more will be enough, depending on the requested precision. 

3 

#10: u := X - X 


# 11 : 

# 12 : 


NEW_ROW(u, X, 0.4472, 10) 

[ 0 . 4472 , -0.4471320317, 0.4467245477, -0.4442912969, 0 . 43009999S2 , -0.3575504623, 0.14S29549S, 
-0.0069S3212209, 6. S1175S641 • 10 , 0, oj 

#13: TABLE(u, X, -1.6, 1.6, 0.001) 
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#15: NEW_ROW(u, X, 0.44725, 10) 

fl6: [0.44725, -0.4474320891, 0.4485269635, -0.4551814749, 0.4984238057, -0.972212844, -1.001238546, 

-1.000002294, -1, -1, -1] 


#17: TAH_FAn(x - X, 0.44725, 10) 



#18: NEW_ROW(u, X, 0.4473, 10) 

#19: [0.4473, -0.4477323983, 0.4503400036, -0.4664759754, 0.5847059419, 15.59091934, 10.40821911, 

6.960229289, 4.672301492, 3.16316677, 2.181452043] 


3 

#20: TAt|_FAI<l(x - X, 0.4473, 10) 



#21: TAr|_FAI<l(x - X, 0.4473, 15) 

#22: NEW_ROW(u, x, 0.4473, 15) 

#23: [0.4473, -0.4477323983, 0.4503400036, -0.4664759754, 0.5847059419, 15.59091934, 10.40321911, 

6.960229289, 4.672301492, 3.16316677, 2.181452043, 1.563843359, 1.207086302, 1.043431923, 
1.002569376, 1.000009843] 
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Expression #8 NEW_ATTR(u , start , end , i nc) returns a list containing all “results” of the Newton- 
Raphson algorithm with initial values xq = start to xq = end with increment inc. So let's try: 


/24: NEWJVTTR(u, -1.5, 1.5, 0.05) 

#25: [-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 0, 0, 0, 0, 0, 


0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , - 1 , - 1 , - 1 , 1 , 


1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 


1 , 1 > 


1 , 


1 , 


1 , 1 ] 


In #9 I define NEW_ATTR_GR(u , start , end , i nc , n) to return a list of points [xq, value of the process 
after n iterations] with start < xq < end. I plot NEW_ATTR_GR(u ,-1.5,1.5,0.05, 10) . (Don't forget to 
activate Approximate before plotting in the Plot Options!) 


#26: NEW_ATTR_GR(u, -1.5, 1.5, 0.05, 15) 



Then we get closer - zooming into the interval [-0.6, -0.4] using an increment of 0.02 - and then into 
the interval [-0.451, -0.446] with inc = 0.00005. 


#27: NEW_ATTR_GRCu, -0.6, -0.4, 0.002, 15) 


- 0.6 


-0.5S -0.56 


-0.54 


-0.52 


-0.5 


-0.4S -0.46 -0.44 -0.42 -0.4 


#28: NEW_ATTR_GR(u, -0.451, -0.446, 0.00005, 15) 


2 -0.4508 -0.4504 -0.45 -0.4496 -0.4492 -0.4488 -0.4484 -0.448 -0.4476 -0.4472 -0.4468 -0.4464 -0.446 


By zooming in we can observe - unexpected - similar patterns. A relationship to “CHAOS” is obvi- 
ous. Selfsimilarity is a significant characteristic of some chaotic evolutions!! 


It is clear that functions with many zeros will show much more interesting patterns. We should ask the 
students to define functions, make investigations and then report about their findings. 

I add one much more striking example: 
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Example 2: 1 want to investigate the behaviour of this algorithm with another function. I am taking 
u = sin 4x, having in mind that there will appear “many” zeros. 


#29: u := SIN(4*x) 

#30: NEW_ATTR_CR(u, -3, 3, 0.05, 15) 

#31: 0.9 ^ X ^ 1.5 

I zoom into the interval [0.9, 1.5] and dis 
cover an unexpected pattern: 


#32: PrecisionDigits := 20 

/33: NotationDigits :i 20 

#34: NEW_ATTR_CR(u, 0 . 9, 1 . 5 , 0 . 005 , 15) 

#35: 1.1 ^ X ^ 1.3 

Next zoom in into [ 1 . 1 , 1 .3] showing another 

pattern. CHAOS IS PERFECT. 8 s 0^9 0.95 1 1.05 1.1 1.1s 1.: 1 . 2 s lA 1.35 1^4 1.45 i.s 1 



tie: NEW_ATTR_GRCu, 1.1, 1.3, 0.001, 15) 

This is an example how I try using the com- 
puter and DERIVE teaching mathematics: 
doing things which I was not able to do ever 
before. 

It is my strong belief that the concepts and the 
ideas are within the capabilities of our stu- 
dents and I hope to transfer another - a new - 
view what mathematics can be. 


[ 1 ] DERIVE - ein mathematischer Assistant, Dr. H.J. Kayser u.a., Landesinstitut fur Schule und Weiterbil- 
dung, 

Soester Verlagskontor 

[ 2 ] Dynamische Systeme und Fraktale, K. H. Becker & M. Dorfler, 1989 , Vieweg, Braunschweig 
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Extending Algebra Concepts 
with Technology {DERIVE} 

Edward Y. Sawada, Mililani, Hawaii 

I recently attended the NCTM conference in Seattle, Washington. I was fortunate to hear Tan Teague 
speak on some of his favourite problems using technology. He illustrated some interesting patterns 
occurring when the coefficients A, B, C vary in the linear equation Ax + By + C = 0. Dan showed 
what happens when A, B and C vary in an arithmetic progression and in a geometric progression. I 
will add to these patterns by showing what happens when A + B = C and A-B = C, B-A = C, and 
V(A^ + B^) = C. Other yatterns will be found by adiustins the coefficients of a Linear Equation . 



Arithmetic Progression 

Linear equation whose coefficients are in an arithmetic progression 
a*x + (a+d)*y + (a+2d) = 0 
Solving for y: 

y = - (a*x + a + 2d)/(a + d) 

Creating a Family of Lines in an Arithmetic Progression: 
VECTOR(VECTOR(- (a-x + a + 2-d)/(a + d), a, -3, 3), d, -3, 3) 

See Fig. 1 


Geometric Progression 

Linear equation whose coefficients are in a geometric progression 

X + r*y + = 0 

Solving for y: 
y = - (x + rVr 

Creating a Family of Lines in a Geometric Progression: 
VECTOR(-(x + r^2)/r, r, -3, 3, 0.25) 

See Fig. 2 



Fig. 1 


Fig. 2 
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Proof that all linear equations having arithmetic progessions for 
their coefficients will intersect at (1, -2). 

We can do the proof with PAP (Paper & Pencil) or we leave the work for DERIVE: 

SOLVE(a*x + (a+d)*y+a+ 2*d= 0Ae*x+ (e+ f)*y+e+ 2*f= 0, [x, y]) 

X = 1 A y = -2 

Proof that all linear equations having coefficients 
Ax + By + (A + B) = 0 will intersect at (-1, -1). 

We can do again the proof with PAP (Paper & Pencil) or we leave the work for DERIVE: 

SOLVECa*x + b*y + (a + b) = 0 A c *x + d*y + (c + d) = 0, [x, y]) 

X = -1 A y = -1 

Other family of equations such as Ax + By + (A - B) = 0 will intersect at (-1, 1) and 
Ax + By + (B - A) = 0 will intersect at (1, -1). These facts can be proven in a similar way. 


Additional Questions from 2006: 

(1 ) Is there any similar family of linear equations which intersect at (1 , 1 )? 

(2) Is there an easier way to prove the conjectures? 

SUBSTCa*x + (a + d)<y + a + 2*d = 0, [x, y], [1, -2]) = true 


SUBST(a*x + b-y + (a + b) = 0, [x, y], [-1, -1]) = true 


ARITHMETIC PROGRESSION 

Examples: 1 , 2, 3, 4, ... 

6, 9,12,15,... 

5, 13,21,29, ... 

4, 14, 24, 34,... 

Examples of equations having arithmetic progression for coefficients: 

1. x + 2y + 3 = 0 2. 6x + 9y+12 = 0 

3. 5x+13y + 21=0 4. 4x+14y + 24 = 0 


Solving any pair of equations for x and j: 

Eq 2 & Eq 4 => (x = 1 , y = -2) 

Eq 1 & Eq 3 ^ (x = 1, y = -2) 

All linear equations having arithmetic progressions for 
their coefficients will intersect at (1, -2). 


Plotting the family of equations on the Tl is possible but it is very time consuming. Notice the function 
for converting the matrix of expressions into a plotable list! 




'yUx) 

/-(x + 3) -2(x + 4) 

I 2 5 

./ -(x + 3) -2 (x + 4) 


-(x + 7) 
4 

-(x + 7) 


-2 


MAIN RflD flUTD FUNC S/ S 









MAIN RAD AUTD FUNC l3iRfl 
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Systems of Equations {Set 1} 

Try solving any two of the following linear 
equations: 

1. x + 2y + 3= 0 

2. 6x + 9y + 12 = 0 

3. 5x+ 13y + 21 -0 

4. 4x + 14y + 24 = 0 

5. -3x + 5y + 13 = 0 

6. 10x-y-12 = 0 

7. X + 3y/2 + 2 = 0 

8. -3x+y + 5 = 0 

9. X - 2y - 5 = 0 

10. 0.1x + 0.8y+ 1.5 = 0 

They all will intersect at (1, -2). 


Systems of Equations {Set 3} 

Try solving any two of the following linear 
equations: 

1. x + 2y-l=0 

2. 6x + 9y - 3 = 0 

3. -X - 2y + 1 = 0 

4. 2x-y + 3 = 0 

5 . 3x + 2y + 1 = 0 

6. 4x + 5y - 1 = 0 

7. X + 3y/2 -1/2 = 0 

8. -3x+y-4 = 0 

9. X - 2y + 3 = 0 

10. 0.1x + 0.8y-0.7 = 0 

All linear equations in the form 
Ajc + By + (A - B) = 0 will intersect 
at(-l, 1). 


Systems of Equations (Set 2} 

Try solving any two of the following linear 
equations: 

1. x + 2y + 3= 0 

2. — 2x + 3y + 1 = 0 

3. 2x- 5y - 3 = 0 

4. -4x - 2y - 6 = 0 

5. 3x-7y-4 = 0 

6. 5x + 3y + 8 = 0 

7. X + 3y/2 + 5/2 = 0 

8. -3x-7y- 10 = 0 

9. x-y = 0 

10. 0.1x + 0.8y + 0.9 = 0 

All linear equations in the form 
Ax + Bj + (A + B) = 0 will intersect 
at (-1,-1). 

Systems of Equations (Set 3} 

Try solving any two of the following linear 
equations: 

1. x-y-2 = 0 

2. 2x + y - 1 = 0 

3. -2x-3y- 1=0 

4. 3x + 2y-l=0 

5. 4x-3y-7 = 0 

6. 7x + 2y - 5 = 0 

7. X + 3y/2 + 1/2 = 0 

8. 9x-7y-16 = 0 

9. 9x + 7y-2 = 0 

10. 0.1x + 0.8y + 0.7 = 0 

All linear equations in the form 
Ax + By + (B - A) = 0 will intersect 
at (1,-1). 


Linear Equation of the form 
Ajc + By + V(A^ + Bb = 0 
will all be tangent to a circle x^+/= 1 

Creating a Family of Lines of this form 
y = - (a X + SQRT(a^ + b+/b 



VECTOR(VECTOR(- (a-x + V(a^2 + b^2))/b, 
a, -3, 3, 1/4), b, -3, 3, 1/4) 
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FAMIL Y OF EQUA TIONS HA VING (a, b) AS ITS SOLUTION 

We have discussed patterns in the linear equation Ajc + Bj + C = 0, and have had various solu- 
tions to the system of familiy members. I would now like to show you how to generate a family of 
equations having the solution (a, b), where you determine the values of a and b. Let the solution dic- 
tate the pattern of the equation. 

If a system of linear equations are to have a solution (a, b), then the linear equations must pass 
through the point (a, b). Since two points determine a line, let (u, v) be the other point the linear equa- 
tion passes through. Using the point-slope form for the equation line, we can generate a pattern of 
family members having (a, b) as its solution. 

E. Sawada's Proof: 

V - b 


#12: 

y - b 

u - a 

#13: 

EXPAM) - b = 

#14: 

a.y + b*u - u*y 

#15: 

(a*y + b*u - 

#16: 

0 = x«(b - v) + 

A = b - V and B = u - a = 
we edit A = a_ and B = b_ 

#17: 

0 = X.(b - (b - 

#13: 

0 = a_*x + b_*y 


(x - a) 
V - b 


i(x - a) .(u - a)j 


which is in our notation: 

Ax + By - (A a + B b) = 0 

2006 Comment: It's much easier to work in vector form: 


^x-a^ 





-0^i^x-a)A + i^y-b)B = 0^ Ax + By- (Aa + Bb) = 0 


or: y-b = [x-a)^B(^y-b) = A(a-x)^ ... 

B 

We intersect two lines of this kind: 

SOLVECp^x + q*y - (a*p + b*q) = 0 A r*x + s*y - (a»r + b«s) = 0 , [x, y]) 

X = a A y = b 

The algebra shown above illustrates that the solution of the system will hQ{x = a, y = b). 


Try solving any two of the following linear 
equations. If (a, b) is to be our solution then 
C = - (Aa + Bb), where A and B are arbitrary 
constants. For this set of problems we will let 
(2, 3) be the solution to the family of equations 


All linear equations in the form 

Ax + By - (2A + 3B) = 0 will intersect at (2, 3). 


Systems of Equations {Set 3} 

1 . ^~^y- 5 = 0 

2. X + 2y - 8 = 0 

3. -X +y - 1 = 0 

4. x-y+l=0 

5. x-2y + 4 = 0 

6. 3x-y- 3 = 0 

7. x + 3y/2- 13/2 = 0 

8. O.lx + 0.8y - 2.6 = 0 
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The following program written in DERIVE will generate n number of equations passing through the 
point (a, b). When solving any two of these family members the solution will be (a, b). 

sys(a, b, n) == VECTOR([RANDOM(10) (-l)^RANDOM(2) (x - a) + 
RANDOM(10H-l)^RANDOM(2)-(y - b) = 0], z, n) 


Function from 1993 does not work 


Output from 


as satisfying as recent versions: 

1993 (V. 2.5x) 


2006 (V. 6.10) 

We experience an unintended sys(i, - 2 , 10 ) 

simplification (sometimes). The 

SYS(1, -2, 10) 

constant is on the right hand side 

9.x - 8. y - 25 = 0 


X = 1 

- or not. Very often we meet two 

- 2.(x - 4.y - 9) = 0 


X = 1 

equal equations following - the 

- 4«x + 5 «y + 14 = 0 


X - y - 3 = 0 

random number generator works 
too slow for the CPU-velocity. 

7-x + S*y +9=0 
6.(x - 1) = 0 


7.x - 5.y = 17 
S.x + y = 6 

1 tried some “tricks” and finally 1 

3.x-y-5=0 


X - y - 3 = 0 

found a way to avoid autosimplifi- 

3.(x - 1) = 0 


3.x + S.y = -13 

cation applying the seldom used 

-S.(x+y+l)=0 


X - y = 3 

quote operator ‘. 

3 . (2.x + y) = 0 


2 . X + 3 . y = -4 

See the program with an addi- 

3.(y + 2) = 0 


. 2.x - 3.y - S = 0 . 


tional parameter k for rising the 
random range. 

sys_6(a, b, k, n, dummy, a_, b_, aux, list, i. 
Prog 

dummy ;= RArOOn(O) 
list := [] 
i := 1 
Loop 

a_ := (RAhDOMCk) + 1) • C-1)^RA®0M(2) 
b_ := (RAHJOPlCk) + 1) • C-l)‘^RAM)OPI(2) 
equ := a_»x + b_<y - (a_*a + b_«b) = 0 
equ ;= LHS(equ) - RHS(equ) 
list := APPE®([[ ' (equ = 0)]], list) 

If i = n 

RETURN list 
i :+ 1 


sys_6Cl. -2, 20, 5) 


10*x - 15*y - 40 = 0 

6*x - IS^y - 36 = 0 

9*x + 17. y + 25 = 0 

X - 10. y - 21 = 0 

17.x - 10. y - 37 = 0 

You can observe the arithmetic se- 
quence for the terms A, B, C in the 
linear equation Ax + By + C = 0. 


lsys_6Cll, -15, 50, 5) 


# 36 : 


19*x + lS*y + 61 = 0 
10*x + 4*y - 50 = 0 
36.x - 49. y - 1131 = 0 
17.x - 2.y - 217 = 0 
. 33. X - 11. y - 533 = 0 . 


#37: 


/33: 


SOLVE 



19 .x + 13*y + 61 = 0 



10-x + 4*y - 50 = 0 



36*x - 49 . y - 1131 = 0 

. [x, y] 


17.x - 2.y - 217 = 0 



. 33 .x - 11. y - 533 = 0 . 
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Last summer I received a letter from Erich containing some sheets of computer pa- 
per with a DERIVE listing. It's headline was “How to program a Maximum - 
Minimum Problem”. The title made me curious and I followed the functions. At 
first I couldn't believe the title, but it is true, indeed. DERIVE is able because of 
the sophisticated - but quite simple - functions to solve Maximum - Minimum 
problems of students' and teachers' daily life. 

Sometimes occur problems, but they are worth to be investigated. Please notice 
how Erich uses ITERATE to work with the solutions of an equation. 


Maximum and Minimum Problems 

Erich Zott, Modling, Austria 


The first example will be treated step by step demonstrating the functions which are finally used to 
“compose the program”. 


Problem 1: 

Given is the circumference f/ of a window (rec- 
tangle + semicircle). Find the dimensions for r 
and h to maximize the area A of the window. 

Model: 

A = + 2rh = Maximum (main condition) 

U = 2r + 2h + rn (side condition) 



/l: 


2 

r *7T 

a := + 2*r*h, s := u - 2*r - 2*h - r*n 


Solve the side condition for variable h: 


#2: ITERATECs, h, 0, -1) = 


u - r* (tt + 2) 


and we substitute into the goal function: 


: 1 im a = 

h^ITERATECs, h, 0, -1) 


r» (2*u - r* (tt + 4)) 


We can use SUBST nowl 

r* (2»u - r* (tt + 4)) 

#4: SUBST (a, h, rTERATE(s, h, 0, -1)) = 

2 

We obtain a function with only one variable a(r) (= F_ONE), 
which then can be differentiated with respect to r: 
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#19: EXTREH(goalj side, r, y, extp) := 


X = EXTR_OHE(ooal , side, x, y) y = EXTR_TWO(goal , side, x, y) 
extp = EXTRVALCgoal , side, x, y) EXTR_DEC(goal , side, x, y) 

d 

F_0HE(goal, side, x, y) — F_OHE(goal, side, x, y) 

dx 


Problem once more: 


# 20 : 


2 

p -n 

s := 1 - 2-p-h, s := ii - 2-p - 2-h - p-n 


II 


#21: EXTT?EB(a, s, p, h) = 


n -I- 4 
2 


h = 


n -I- 4 


II 


extp = 

2- (n -1-4) 

p- (2- II - p- (n ■+■ 4)) 


RaxiHii 


II - p - (n -I- 4) 


Let's change the order! (Which is solving the side condition for h and proceeding ....) 


#22: EXTPEII(a, s, h, r) = 


h = 


n -I- 4 
2 


n -I- 4 


extp = 


2- (n -1-4) 

(ii - 2- h)' (2- h- (n -I- 4) + n-ii) 
2 

2-(n + 2) 


HaxiHii 


4- (ii - h- (n -I- 4)) 


(n + 2) 


All the sketches are scanned from the original DNL#12 from 1993. 


Problem 2: 


m 




i f- 



9<foo 


yyi 




2tr) 


Mir)inwnt [ 


Modell: 


z = (x + 4)(_y + 2) -x_y = Minimum 
xy = 9800 


#1 : LOAD (D : \DfD\DNL\DNL93\nTHl 2\naxn i n_ut . mth) 

#2: [water := (x + 4)»(y + 2) - x»y, ar := x»y - 9800] 

#3: EXTREn (water , ar, x, y) 

X = 140 y = 70 

extr = 568 ninimum 


#4: 


2 

2.(x + 4.x + 19600) 


X 


2.(x - 19600) 

2 

X 


X = 140 and 7 70 make a minimum shaded area 

of 568 square units. 
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Problem 3: 



The triangles are equilateral. 

Modell: 


area = ab + 2b^ 


2a' V3 . 

H = Maximum! 

4 


20 = 4a + 6b 


#1 : LOAD(D : \DfD\DNL\DNL93\nTH12\naxnin_ut . mth) 

2 

2 2. a .^3 

#2: area := a«b + 2>b + 

4 


#3: 

#4: 


#5: 


circf := 20 - 4*a - 6*b 
EXTREn(area, circf, a, b) 
450.^3 200 


227 227 

5600 1250.^3 

extr = 

227 227 



300.^3 

227 


ninimum 


2 

a *(9.V3 + 4) - 100. a + 400 
18 

a = 2.552523627 
extr = 15.13187881 


a. (9.^3 + 4) - 50 
9 

b = 1.631650915 
ninimum 


# 6 : 


2 

a .(9.^3 + 4) - 100. a + 400 a. (9.^3 + 4) - 50 


L 18 9 J 

The interesting answer is: the area of 15.13 is a 
Minimum! Check it by plotting T ONE, (bottom 
left) 


Problem 4: 



Vol = u^n + 


(4- v)7T 
3 


= Maximum! 


Av = \6-u^ 


#i: a 


OAD(D : \DfD\DNL\DNL93\l'ITH12\naxnin_ut . mth 


# 2 : 


#4: 


#6 


2 

2 U *77 2 

vol := u .77*v + *(4 - v), par := 4.v + u - 16 


#3: EXTREf1(vol, par, u, v) 

u = 0 V = 4 

extr = 0 ninimum 


2 2 2 
77*U *(24 - u ) 2*77.u.(12 - u ) 


6 3 

#5: EXTREn(vol, par, v, u) 

V = 1 u = 2.^3 

extr = 24.77 Haximum 

8*77* (v + 2). (4 - v) 16*77.(1 - v) 


Problem 5 is quite interesting because it worked properly with DERIVE 2.xx and DERIVE 
3. XX (DOS DERIVE) as shown below. I reactivated my DOS-DERIVE and checked it - it 
really works. 
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This is a part of the DERIVE 6 output. 

y = -3 


X = 0 V X = 0 


extr = 0 


'/(s - y )*Cy + 3) 


2 2 

IF(?, naximum, IF(DER_2hD(x< (y + 3), x + y - 9, y, x, 1) ; 

2 y.(y + 3) 

V(9 - y ) 


VC9 - y ) 


The reason for the problem lies in another internal programming of the ITERATES-function. 
In EXTR_ONE we solve an equation which can have more than one solution - usually in 
school problems might appear quadratics. As you can see I changed EXTR_ONE applying 
SOLUTIONS with an additional parameter k_ with defaults to 1 (1®‘ solution) or must be en- 
tered as 2 (2"'' solution). Then I adapted all other functions accordingly 
(see file MaxMin_ut_impr.mth). 

S — F_ 


# 2 : 


EXTR_ONECa, s, r, h, k_ := 1) :i j^SOLUTIONSj^— F_ONE(a, s, r, h) 


-]] 


Then I tried once more using the improved version: 

#1 : L0iAD(D:\DfD\DNL\DNL93\HTm2\Haxl1iii_iit_iHpr.Hth) 


# 2 : 


trap_ap 


o + 2-y + 3)x 2 2 

, pyth ;= 9 - X - y 


#3: EXTRH(tpap_ap, pyth, x, y) 


3-^3 


X = 


extr 


2 

9-^3 


y = - 


RaxiBii 


x.(J(9 - X ) + 3) ^(9 - X ) - 


+ 3 


2 

^(9 - X ) 
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#4: EXTREIlCtpap_ap, pyth, x, y, 2) = 


3-^3 


2 

9-^3 


y = - 


extr = - 


2 2 X 

x-ao - X ) + 3) 7(9 - X ) 


7(9 - X ) 


+ 3 


#5: EXTREll(tpap_ap, pyth, y, x) 


y z -3 


extp = 0 IF(?, HaxiHiH, IF(DER_2ND(x- (y + 


7(9 - y )-(y + 3) 


#6: EXTREll(tpap_ap, pyth, y, x, 2) - 


3-73 


2 

27-73 


extr = “ 


7(9 - y )-(y + 3) 7(9 - y)- 


Haxiaii 


2 y(y + 3) 


7(9 - y ) 


The negative x-value is caused by the ITERATION-function - there are two solutions for the 
underlying equation and ITERATION returns now the negative one. Don't know why. Prob- 
lems occur in all cases with more than one solution. 

We see that the first three attempts are leading to solutions which are outside of the range 
because of the negative y-value. The fourth try is successful. 

Much later, in 2002 in wrote a book “Programmieren in DERIVE”^^'. I came back to Erich's 
basic ideas and tried to program the MaxMin-procedure with DERIVE 6: 


2 2 

extreiiCx'O + y)i3f +y = 9, y, x) 


y X Optniiuii 

1 1-4^ 27--/3 

2 2 4 

1 1-4^ 27--/^ 

2 2 4 

. -I 0 0 


Art 

1oka1es llniniiuii 

1oka1es naxniiuii 
kenn EW oder nicht def. 


Josef Bohm, Programmieren in DERIVE, bk-teachware Schriftenreihe SR-32, bk-teachware 2002 
(For all of you who are owning this “masterpiece^ is an improved version of program extrem included 
among the DNL 12-files.) 
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In 1993 Keith Eames gave step by step instructions how to work with Derive 2.xx. 
So he wrote: 


Step 1: <A> r'^S <C> <S> <ENTER> <ENTER> <ENTER> 

I'll try to keep his style In best possible form but adapt It for the recent Derive ver- 
sions. Josef 

Use of DERIVE in Proofs by Induction 

Keith Eames, Romford, UK 


Assumed Knowledge: Derive commands, for example SUM, VECTOR, TABLE, de- 

claring functions, I-notation. Introduction to method of Proof by 
Induction. 

The following two examples will illustrate how Derive can help to reinforce your understand- 
ing of the Method of Proof by Induction. 


Example 1 

Prove by induction f +2^ +3^ +... + n^ 


n^{n + ^f 
4 


Step 1 : <F2 (= Author Expression)> <ENTER> 

<Alt + C> <Sum> <Ok> or 
shortcut <Ctrl + Shift + S> <Ok> 

This will give the series being looked at in summation form. 



Step 2: <F2> S(n):=#2 <ENTER> 

or S(n):=sum(r'^3,r,1,n) <ENTER> 

This will declare a function S which will contain the answers to the sum of for given values 
of n. 


Step 3: <F2>T(n):=n''2 (n+1)'"2/4 <ENTER> 

This will declare a function T which will contain the answers to 1^+2^+3^ + ... +n^ obtainable 
from the quoted formula. 


Step 4: <F2>TABLE([S(n),T(n)],n,1,10) <ENTER and Simplify> 

This will compare the two sets of answers for each value of n in S{n) and T(n). This should 
lead to a comment that the result seems to be true for n = 1 to 10. 

Now assume the statement 

r +2+3 + ... + n = — ^ — 


4 


is true for n. 
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steps: <F2>S(n) + (n+1)^3 = <ENTER> 

By adding the (n+1)th term to the original series we are required to show that 

S(n+1) = S(n) + (n+1)3. 

Step 6: <F2>S(n+1) = <ENTER> 

On looking at the answers to the line numbers for step 5 and step 6 (#7 & #8), the answers 
should be identical. User Der/Ve’s FACTOR command on each answer to verify that the two 
statements are identical (or subtract one answer from the other one to obtain 0). 

Students should be able to recognize the identity without using FACTOR, of course. 

n I 

S(n) := Z r 

r=l 

2 2 
n -Cn + 1) 

#4: TOO := 

4 

#5: TABLECtSOO, TOO], n, 1, 10) 

11 1 
2 9 9 

i 56 56 

4 100 100 

5 225 225 

# 6 : 

6 441 441 

7 7&4 7&4 

a 1296 1296 

9 2025 2025 

. 10 5025 5025 



#10: FACTOR 


4 


4 
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Example 2 

Prove by induction: 7" + 2 is always divisible by 3 for n > 1 . 

Step 1 : <F2 > + 2 <ENTER> 

This will give the numbers being looked at for various values of n. 

Step 2: <F2> T(n):=#1 <ENTER> 

or <F2> T(n):=7''n+2 <ENTER> 

This will declare a function T which will contain the answers to the values of 7” + 2 for de- 
fined values of n. 

Step 3: <F2> VECTOR(T(n),n,1,10) <ENTER and Simplify> ^ giving expression #4 
<F2> #4/3 <ENTER and Simplify> ^ resulting in integers only 
or <F2> mod(#4,3) <ENTER and Simplify> ^ resulting in a list of zero reminders 

This will show whether the conjecture seems to be true for n = 1 to 10. Now assume that the 
statement ist true for n. 

You are required to show that T(n+1) is also divisible by 3. 

Step 4: <F2>T(n+1) <ENTER and Simplify> 

By looking at the expression we have to show that it is divisible by 3. 

Step 5: <F2> 7T(n) <ENTER and Simplify> 

Comparing this answer with T(n+1) we can show that T(n+1) = 7T{n) - 12 which is divisible 
by 3 under the assumption that T{n) ist divisivie by 3. 

The print out of this example is given in the following, with statements added: 

(We can prove this without any CAS: rewrite T(n+1 ) = 7'^*'' + 2 as 7 - 7" + 2 = 6 • 7" +7" + 2 
and this is obviously divisible by 3, because (7" + 2) and 6 are both divisible by 3.) 

n 

fl: 7+2 

n 

#2: TOO := 7 +2 

#3: VECTORaOO. II. 1. 10) 

#4: [9, 51, 345, 2403, 16809, 117651, 823545, 5764803, 40353609, 282475251] 

[9, 51, 345, 2403, 16809, 117651, 823545, 5764803, 40353609, 282475251] 

#5: 

3 

#6: [3, 17, 115, SOI, 5603, 39217, 274515, 1921601, 13451203, 94158417] 

#7: nOD([9, 51, 345, 2403, 16809, 117651, 823545, 5764803, 40353609, 282475251], 3) 

# 8 : [ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ] 
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True for n = 1 to 10 

Assume T(n) is divisible by 3 

Prove thatT(n+1) is also divisible by 3 

#9: Kn + 1) 

II + 1 

#10: 7 +2 

Compare v\(ith T(n) = 7^n + 2, multiplying T(n) by 7: 

II + 1 

#11: 7-TOi) = 7 +14 

II + 1 II + 1 

#12: Kn + 1) = 7-Kii) - 12 = (7 +2 = 7 +2) 

T(n) and 12 are both divisible by 3 then T(n+1) is also divisible by 3. 
qued. 


It is very easy to transfer proofs like these onto the TI-92 or the Voyage 200: 


["?2fco.T3nd1u?™fE>cJ^utefFf^d...1 1 

f"^1mS..fc?fcl0tTer1p.?.I0fcleyn UpI 1 

f"^1fllSp.fc?fcl0tTe4pr?.I0lcie;‘n Upf 1 

C'-K r-^3 , r , 1 , n > -^s C n > 
C:n-^2+(:n+l>^2Y4-^t(:n> 
CiseqCsCk^, k, 1, 1G> 
CiseqaCk^, k, 1, 1G> 
CisCn^+Cn+l^^S 
CisCn+n 

CisCn^ + Cn+n-^S-sCn+O 
C : factor tstn> + tn+n-^3> 
C: f actor tstn+n > 

■ X Done 

r=l 

n^ (n + 1)^ 

■ ^ t(n) Done 

■ seq(s( k) , k , 1 , IG) 

a 9 Z6 IGG 225 441 784 1298 ^ 

■ = onf + f L-h L- 1 ini 

■ seqCUk), k, 1 , IG) 

a 9 38 IGG 225 441 784 1298 ^ 

. -fri 1 fr 1 iP (n+l)^ (n2 + 4 r, + 4) 

■s(.nj+t_n+i_j ^ 

Cr, + 2)2(n2 + 2r,+ l) 

■ s(.n + ij ^ 

sea<t<k>.k.l.lG> 

s<n>+<n+l>^3-s<n+l>l 

MAIN RAD AUTD FUNC 

MAIN RAD AUTU FUNC 5/30 

MAIN RAD AUTU FUNC 5/30 




' sCn + 1) 


■ sCn) +(n + 1)^ - sCn + 1) 

■ factor(s(n) +(n + 1)^) 


G 

(n+ l)^(n + 2)^ 
4 


■ factor(s(n +1)) 

s<n>+<n+l>^3-s<n+l>l 

MAIN RflD flUTD FUNC S/SO 


^■(n + 2)^ 


FT? FH YTe Y" 

\j ICopiPiand |U i eto |Execuie |F i nd... | 


C:7-^n+2H(:nV 
CiseqaCk^, k, 1, 1G> 
CiseqCniodacb^S^, k, 1, 1G> 
CitCn+n 
C -. 6 * 7 ^ n + t (. n >\ 


] 


MAIN RflD flUTD FUNC 


Use a CAS to help proving the following statements by the Method of Induction: 

(a) 2^"“^ + 3 is divisible by 7 (b) Show that 6 is a factor of + 5n 


(c) ^rr! = (n + 1)!-1 
1 n 


r=1 


(e) E 
(g) 


l■^7^ 4r^ - 1 2n + 1 


1 1 

+ - 


■ + ... + ■ 


1 


1 


(d) 3^"*^ - 8n - 9 is divisible by 64 

(f) ^4r-3^ =3"^'(2n-1) + 3 

r=^ 

1 


1-2-3 2-3-4 n(n + 1)(n + 2) 4 2(n + 1)(n + 2) 

(h) 1 +5 + 9 + ... + (4n-3) = n(2n-1) 

1 n ^ ^2 _ n{n + ^){2n + ^) 

/(=1 


(i) 


1 1 
■ + + ... + - 


1-2 2-3 n(n + 1) n + 1 

3 factor of 5‘ 
(m) 1 + 2 + 2^ + ... +2"-^ =2"-1 


(k) 


(I) Show that 1 9 is a factor of 5^”"' • 2""" + 3""" • 2^""' 

■>2 _L _L ryn-^ — on A /«\ ^f\n+^ 


(o) 1-5 + 2-5" +3-5^ +... + A7-5" = 


(n) 10"^' -9n- 10 is divisible by 81 

5 + (4n-1)-5”^' 

16 
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I am very glad to welcome the first Asian contribution in the DNL. We met Mr Chuan at the TMT 
Conference in Birmingham. He surprised the participants of the UK-DUG Meeting with his enthu- 
siastic talk about using DERIVE and a spreadsheet program. He promised to submit papers for the 
DNL. Thanks to Taiwan. 

m ^ ^ ^ ^ ^ ^ 

NATIONAL TSING HUA UNIVERSITY 

DEPARTMENT OF MATHEMATICS 

HSINCHU, TAIWAN 30043, REPUBLIC OF CHINA 
TEL: (035) 713704 

Playing Cards Shuffling with DERIVE 

Jen-Chung Chuan, Hsinchu, Taiwan 

When a deck of 52 playing cards are perfectly shuffled, there associates a permutation of the form 

1 2 3 4 ... 49 50 51 52" 

1 27 2 28 ... 25 51 26 52^ 

describing how the position has been altered. The iterated k fold product 



k 


therefore represents the position of the cards after identical shuffling has been applied ^-times. The 
following three lines of inputs in DERIVE computes /, r^, ... until the identity permutation has 

been reached. 

u := VECTOR Cn, n, 52), k := IFlflOOCn, 2), + 26, M 

SHUFFLECv) := VECTOR (v , n, 52) 

#2: k 

ITERATES (SHUFFLE Cv), v, u) 

Does it surprise you that the cards are back to their original position after shuffling only 8 times? 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

1 

27 

2 

28 

3 

29 

4 

30 

5 

31 

6 

32 

7 

33 

8 

34 

9 

1 

14 

27 

40 

2 

15 

28 

41 

3 

16 

29 

42 

4 

17 

30 

43 

5 

1 

33 

14 

46 

27 

8 

40 

21 

2 

34 

15 

47 

28 

9 

41 

22 

3 

# 4 : 1 

17 

33 

49 

14 

30 

46 

11 

27 

43 

8 

24 

40 

5 

21 

37 

2 

1 

9 

17 

25 

33 

41 

49 

6 

14 

22 

30 

38 

46 

3 

11 

19 

27 

1 

5 

9 

13 

17 

21 

25 

29 

33 

37 

41 

45 

49 

2 

6 

10 

14 

1 

3 

5 

7 

9 

11 

13 

15 

17 

19 

21 

23 

25 

27 

29 

31 

33 

. 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 




Part of the DERIVE output 




